ROTATING ELASTIC BODIES IN EINSTEIN GRAVITY 



LARS ANDERSSONt, ROBERT BEIG*, AND BERND G. SCHMIDT 

Abstract. We prove that, given a stress-free, axially symmetric elastic 
body, there exists, for sufficiently small values of the gravitational constant 
and of the angular frequency, a unique stationary axisymmetric solution to 
the Einstein equations coupled to the equations of relativistic elasticity with 
the body performing rigid rotations around the symmetry axis at the given 
angular frequency. 
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In the paper [T], we constructed for the first time static, self-gravitating elas- 
O^. tic bodies in general relativity with no symmetries. Here we build on the ideas 

}L( \ and techniques introduced in that paper to construct solutions to the Einstein 

equations describing steady states of self-gravitating matter in rigid rotation. 
The matter model we use is, as in [T], that of a perfectly elastic solid. We 
make the minimal symmetry assumptions necessary for a steady state in rigid 
rotation, namely we assume that the reference body has an axis of symmetry. 
I Further, we assume that the elastic material is isotropic. This condition, which 

^\ ' was not needed in the static case, is necessary for our construction in the case 

of a rotating body. 

The only class of solutions to the stationary Einstein equations with rotating 
matter previously known are the rotating perfect fluid solutions constructed by 
I Heilig for a certain class of equations of state. In the Newtonian theory, 

existence of steady states of self-gravitating perfect fluids in rigid rotation was 
. established by Lichtenstein, see [8] for a modern presentation, and by Beig and 

rS I Schmidt [5] for the case of elastic matter. All these solutions, including the 

■ ones constructed in the present paper are in addition to being stationary, also 

axisymmetric. 

In the Newtonian theory, two families of non-axisymmetric rotating fluid 
configurations in explicit form are known, see [15] and references therein. These 
families of solutions are the Dedekind ellipsoids, and the Jacobi ellipsoids, which 
in the language of general relativity have helical, but no stationary or axial 
symmetry. 

One expects asymptotically flat rotating solutions of the Einstein equations, 
which are not axially symmetric to be radiating, and hence non-stationary. 
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However, if one relaxes the condition that asymptotic flatness holds in the 
usual sense, it may be possible to construct helically symmetric solutions of the 
Einstein equations which are not axially symmetric. See [2] for a study helically 
symmetric solutions in the special relativistic case. An argument to the effect 
that axisymmetry necessarily holds for a rotating fluid in general relativity was 
given by Lindblom [T?], assuming that the fluid is viscous. 

Equilibrium states of fluids or collisionless matter play an important role in 
astrophysics, providing the basic models of stars and galaxies. Depending on 
the equation of state, or in the case of collisionless matter, on the properties 
of the distribution function, a steady state may describe a compact body, or 
a configuration where the matter density is nowhere vanishing. Typically, the 
objects of interest are compact. 

In addition to fluids and collisionless matter, elastic bodies are of considerable 
interest in astrophysics in view of the fact that there are strong theoretical 
reasons for supposing that neutron stars have a solid crust, modelled by elastic 
matter, cf. e.g., [iniEllT]. The solutions which have been constructed in the just 
mentioned papers are all spherically symmetric, although perturbation analyses 
have been carried out, allowing for axial perturbations breaking the spherical 
symmetry 

1.1. Rotating bodies in elasticity. Elastic matter is, as discussed in section 
12. II below, described by a map from spacetime to a three dimensional mani- 
fold, called material manifold or body, whose points label the particles making 
up the elastic continuum, and which is taken to be a connected, bounded do- 
main in flat R^. 

In considering a rotating steady state, it is important to distinguish between 
the microscopic and macroscopic degrees of freedom. The microscopic degrees 
of freedom of the elastic matter are described by the configuration f^, while 
the macroscopic aspects are described by the stress energy tensor generated 
by the matter, and the metric of the spacetime containing the body. For a 
rotating body in equilibrium, it is the case that the stress energy tensor, as 
well as the spacetime metric are stationary, i.e. invariant under the flow of a 
Killing vector ^^d^ = dt, called the stationary Killing vector, while the matter 
particles, described by the configuration map, are in motion relative to dt- 

As mentioned above, the rotating bodies we construct are axially symmetric. 
In particular, the spacetime containing the body admits a Killing field r/'^S^ = 
dfj), called the axial Killing vector, which commutes with In addition, there 
is a constant the angular frequency of rotation, such that the matter particles 
move along the helical orbits of the Killing vector ^^-|-J7r/^, i.e. the configuration 
is constant along the fiow of the helical field. 

It is nevertheless the case, assuming axisymmetry of the configuration, and 
that the elastic material is isotropic and frame indifferent, cf. section 12. H 
that all spacetime tensors naturally derived from the configuration are both 
axisymmetric and stationary. This applies for example to the matter flow vector 
induced by the configuration, the stress tensor or, in fact, to the full stress 
energy tensor. 
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We now briefly describe the method used in this paper. Consider a Cauchy 
surface M transverse to the stationary Kilhng field. The equations for the 
gravitational field variables are derived by imposing the condition of stationarity 
and restricting the Einstein equations, reduced in harmonic gauge, to M. No 
axisymmetry condition is imposed at this stage. 

The Einstein equations imply, via the Bianchi identity, a set of equations for 
the matter variables. These equations are derived by considering a configuration 
which is comoving with respect to the helical flow, and restricting to M. Here 
the axial vector field ry^ is assumed to be specified in advance. There are, a 
priori, four matter equations for the three unknowns f"^. We deal with this 
problem by simply dropping one of the four equations. It turns out, however, 
that this supplementary equation follows from the others when r/^ is Killing, as 
is the case for a solution to the system derived by the above procedure. 

The resulting coupled system of equations, assuming standard constitutive 
conditions for the elastic material, is elliptic for suffiently small values of Q,. 
The system depends on the parameters G and $7. We look for solutions to this 
system for small, nonzero, values of G, near the background solution given 
by taking the spacetime to be Minkowski, the configuration to be stressfree and 
the Newton constant G and $7 to be both zero. 

The boundary between the matter region and the vacuum region depends on 
the unknown configuration. To deal with this problem we write the equations in 
material form in a way analogous to [Ij and apply the implicit function theorem. 
This can not be done directly due to the failure of the linearized operator to 
be invertible. This is a standard problem for elasticity with natural boundary 
conditions, i.e. vanishing normal stress at the boundary. Following [l], what 
we actually solve is a projected version of the field equations, such that the 
implicit function theorem does apply. One must then show, as in fact turns out 
to be the case, that the solution to this projected system is actually a solution 
to the full system. 

So far the vector field rj^^ was essentially arbitrary except for the condition 
that it commute with In order to ensure that r]^ is a Killing vector, we now 
assume that the material manifold is axisymmetric as a subset of Euclidean 
and that r]^ is the pull back of the axial vector field on the body under the 
trivial configuration. It then follows by uniqueness that the vector field i]^ is a 
Killing vector. 

It now remains to show that the solution found by the implicit function ar- 
gument satisfies the Einstein equations. In particular, the gauge conditions 
must be satisfied and the elastic equation must be valid in its original form. 
This condition is equivalent to the condition that a certain linear system of 
equations coming from the Bianchi identities has only the trivial solution. In 
fact, the linear system under discussion is homogenous precisely because the 
Killing nature of r]^^ guarantees that the above mentioned supplementary equa- 
tion is satisfied, provided that the main elastic equation is valid. The rest of 
the argument follows essentially the pattern of [1]. 

1.2. Outline of the paper. In section [2T| we give some background on rel- 
ativistic elasticity. Section 12.31 introduces stationary metrics and defines the 
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field variables hik,U,ipi used to parametrize the spacetime metric. In contrast 
to the static case, there is a further component ij^i of the gravitational field, 
corresponding to the failure of to be hypersurface orthogonal in general. 
Next, in section 12.41 the stationary Einstein equations are written in terms of 
the field variables just introduced, this corresponds effectively to performing a 
Kaluza-Klein reduction. 

The field equations imply a set of integrability conditions, which are derived 
in section 12. 5[ One of these identities is the elasticity equation, which is later 
used as one of the set of equations to be solved using the implicit function 
theorem. The rotation of the elastic body is introduced in section 12. 6i This 
is done by choosing a spacelike vector field rj^^d^ which commutes with the 
stationary Killing vector and assuming that f^,^ (^^ + r^r/^) be zero. It will 
later turn out that jy'^ is actually the axial Killing vector. 

In section [2^1 the stress tensor is expressed in terms of the geometric variables 
and the configuration f^. In particular, this allows us to write the components 
of the stress energy tensor in terms of the field variables, and obtain (I2.37p . the 
basic PDE system for hik,U,ilJi, f^- As some of the equations are not elliptic 
we use, as in the static case, harmonic coordinates to extract an elliptic system. 
In the stationary case, it is necessary to make explicit use also of the condition 
that the time function be harmonic, cf. section 12. 8[ Finally, we have reduced 
the field equations to an elliptic free boundary value problem in space. In order 
to avoid dealing directly with the free boundary aspect of the problem, we 
move all equations to the body. This is done in section 12.91 following closely 
the procedure in [T]. The final details needed to completely specify the PDE 
problem to be solved are introduced in section 12.101 There we introduce the 
relaxed state and a fiat metric on the body. We assume that the shape of the 
body is axi-symmetric and use the relaxed configuration to define the vector 
field ??* in space. We also introduce the assumption that the elastic material is 
isotropic. 

As in the static case considered in [1], we must consider a projected system in 
order to be able to apply the implicit function theorem. The analytical aspects 
of this problem are considered in section [3l The solution to the projected system 
is then shown in section H] to be a solution to the full set of field equations in 
the body frame, and to be axisymmetric. In section we move the projected 
equations and solutions back to space. The vector field r/* is proved to be a 
Killing field in section 14.21 

In section 14. 3^ we derive some divergence identities play an essential role in 
the equilibration argument. This leads up to the main theorem 14. 6( which is 
stated and proved in section 14. 4[ In particular, we prove that the harmonic co- 
ordinate conditions are satisfied for the solution of the reduced system that has 
been constructed, and hence that we have solved the full set of field equations. 

The spacetimes constructed in Theorem 14.61 have an isometry group M x S"^ 
generated by the commuting Killing fields S^^,r]^. For spacetimes with two 
commuting Killing fields, it was first proved by Papapetrou [T7] in vacuum 
and by Kundt and Triimper jl3] for fiuids, that orthogonal transitivity holds. 
Recall that orthogonal transitivity is the condition that the distribution of 
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2-surface elements perpendicular to the generators of the symmetry group is 
surface forming. This condition is used for constructing Weyl-type coordinates 
which play a dominant role in attempts in the exact solution literature to find 
rotating body solutions, see [18] and references therein. 

In section [^31 we establish that for the spacetimes constructed in theorem l4.61 
the distribution defined by the 2-surface elements orthogonal to the group orbit 
for the action of the stationary and axisymmetric Killing vector is integrable. In 
the case of a smooth spacetime, where the Frobenius theorem applies directly, 
this fact implies that orthogonal transitivity holds, i.e., that there are 2-surfaces 
perpendicular to the generators of the symmetry group. In the present case, 
however, this step needs further analysis which we defer to a later paper. 



2. The field equations of a rotating, self-gravitating elastic 

BODY 

2.1. Relativistic elasticity. Let {A4,g^u) be a 3+1 dimensional spacetime. 
The body ;B is a 3-manifold, possibly with boundary. We shall consider the case 
when ;B is a bounded domain in the extended body The body domain B is 
assumed to have a smooth boundary. In this paper, we shall only consider the 
case where B is connected. The fields considered in elasticity are configurations 
f : M ^ B and deformations <j) : B ^ A4, with the property that f o (p = idg. 

Let t be a time function on M and introduce a 3+1 split = R x M. 
We consider coordinates (x^) = (t, x*) on A^, where are coordinates on the 
space manifold M. On we use coordinates X"^. The body B is endowed 
with metric 5 13 and a compatible volume form Vabc- We assume that in a 
suitable Cartesian coordinate system 5b has components 5ab where 8ab is the 
Kronecker delta, and V123 = 1. 

The configuration f : M ^ B ishy assumption a submersion. The derivative 
of / is assumed to have a timelike kernel, i.e. there is a unit timelike vector 
field on f~^{B) with u^u'^ = —1, such that 

The field is the velocity field of the matter and describes the trajectories of 
the material particles. 

Let A = A(/, df, g) be the energy density for the elastic material in its own 
rest frame. The Lagrange density for the self-gravitating elastic body now takes 
the form 

^-^-^^^ 

The Einstein equations resulting from the variation of the action with respect 
to g^'^ take the form 

where Gf^^ is the Einstein tensor of gf^^ and T^i, is the stress energy tensor of 
the material, given by 

OA 
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On the other hand, the canonical stress energy tensor is given by 

~ QfA J 'A* • 

General covariance implies, by the Rosenfeld-Belinfante theorem, that 

T — — T 

see [121 section 7]. 

Given a configuration /^(x^), define 7^-^ = f'^^^f^^vQ^^ and let ^ab be the 
inverse of 7"^^. General covariance implies A is of the form A = A(/^,7'^^), cf. 
[121 section 7] , see also [H section 4] . A stored energy function of this form is said 
to satisfy material frame indifference. If in addition, as we shall later assume, 
A depends only on the principal invariants Aj, i = 1, 2, 3 of 7'^'^, defined as the 
elementary symmetric polynomials in the eigenvalues of ■^^b = l"^^ {^b)cBj 
then the material is called isotropic. 

Define 

dk 

Then we have 

T^y = kU^Uy + S^y , (2.1) 

where S^j_y = SAsf^^^if^^v In particular S^yU^ = 0. The relativistic number 
density Ug is defined by 

ng = —Vabcva'B'C'7 7 7 



We have Ug = (det7'^^)^/^ and hence 

dua 1 

= r^^AB . (2.2) 

Define the stored energy function e by 

A = nge, (2.3) 
and the elastic stress tensor tab by 

de 

^AB = . 
With these definitions. Sab takes the form 

Sab = ngTAB , 

and we can write 

Tfj,y = Ugeu^Uu + ngTABf^,^lf^ ,u ■ 
See [19] for a more explicit expression of T^y in terms of the invariants (Aj). 

If material frame indifference holds, then if A is viewed as a functional of 
f^,gfj,u, we have that for any spacetime diffeomorphism a, 

A[foa,a*g]=A[f,g]oa, 

and hence all spacetime quantities constructed from A^f^^^g^^, are covariant 

ig, u^" and TABf^,iij 



under a, including n„, and tabI^ nf^ v In particular, this holds also for 



T 
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Let S be an isometry of {B, Sb)- The matrix {'E^j)^b is related to j^b by an 
orthogonal similarity transform and hence has the same invariants Aj as j^b- 
Hence, for an isotropic material, 

A[So/,5] = A[/,5]. 

2.2. Material and spacetime isometries. We now introduce the notion of 
symmetry of a configuration which will play an important role in this paper. 
Suppose the spacetime {M,g) has an isometry a. Then a defines a material 
symmetry of if there is an isometry S of {B,Sb) such that 

So/ = /oa. 

Thus, in particular, if the configuration is comoving with an isometry, i.e., if 
is proportional to a Killing vector then the configuration has the flow ag 
of ^'^ as a material isometry. with S given by the identity map on B, in which 
case it follows that f^oOf^ = f"^. However, this does not hold for a general one- 
paramctcr family of material isometrics. It follows from the last two statements 
in the previous subsection that, for an isotropic material, a spacetime isometry 
a which also defines a material isometry leaves the Lagrangian invariant, i.e. 
■^[/'fl'] — ^[f^g] ° ^'^'^ thus T^jy is also invariant under a, i.e. a*T = T. 

The following is an example which is relevant for the situation in this paper. 
Suppose we have two timelike Killing vectors and ^'i^. In the situation 
considered in this paper, the interesting case is where E,^ is the stationary Killing 
field, while ^'^ = ^^ + nri^^ is the helical Killing field. Then one may consider the 
case where the configuration is comoving with respect to E,'^ while the flow ag of 
defines a configuration symmetry in the sense that there is a one-parameter 
family of isometries of the body such that 

In this case. Eg are rotations of the body. We sec from the above that it is 
possible for the configuration to explicitly depend on the Killing time i, defined 
with respect to although T^j_^ is independent of t. 

2.3. Stationary metrics. We now assume {M,g) is stationary, i.e. there is 
a timelike Killing field ^^d^i = dt- Further we assume the space manifold M 
is diffeomorphic to M^. It will sometimes be convenient to denote this space 
by M.g. Define a function U = ^ In^'^^^ and a one- form = i/jidx^ such that 
e~'^^ ^^dx^ = dt + ip. Then g can be written in the form 

g^^ dx^'dx'' = -e^^{dt + ^ndx'f + e'^^ hijdx'dx^ , (2.4) 

where hijdx^dx^ is a metric on the level sets of t, and U,xpi,hij are time inde- 
pendent. The inverse metric takes the form 

g^-d^d, = -e-2^af + e'^'^h'^idi - m){dj - ^jdt) , (2.5) 

where h^^hjh = S^k. The spacetime volume element is given by 

,/^=e-'^^Vh. (2.6) 

The assumption that ^'^d^ = dt is & Killing vector implies 

Bgt = -e^^D^'^i Ugx' = -e^^h^^T]^ , (2.7) 
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where Di and the Christoffel symbols refer to hij. 

2.4. Kaluza-Klein reduction. Let ojij = d^iipjy The scalar curvature Rg for 
a metric of the form ()2.4p is given by 

RgV^ = Vh{Rh + 2AhU - 2\DU\l + e^^\Lo\l) . (2.8) 

Here \DU\l = DiJJD'^U and similarly for \uj\\. Define H^^ by 

The reduced number density n is defined with respect to H^^ , 

= ^VabcVa'B'C'H^^' H''''' H^^' . 

Then we have 

dfi 1 

Qj^ = 2^^^^ ' (2.9) 

ng = e^^n, (2.10) 
so that with the form ()2.4p for g, we have 

A^/^ = nee^Vh. (2.11) 

Taking into account the fact that the term 2AhU in the scalar curvature ex- 
pression (j2.8p contributes a total divergence to the action and can be dropped, 
we may now write the action in the reduced form 

C = (Rh - 2\DU\l + e'^Hl) +pe^Vh, 

where p = ne. Let Gij = Rij — \Rhij be the Einstein tensor of hij and define 
= -^[{D,U){D,U) - hi, iDkU)iD'U)] , (2.12) 

and 

^ij = ^ e'^^ih,, UMUJ''' - UikUJj''] . (2.13) 
The reduced field equations now take the form 

AhU = 47rGe^ (^p + ^) Xf-^B) " e^^'^fc/u;"' , (2.14a) 
D^(e^%) = -87rGe^|^X/-i(B) , (2.14b) 
Gij = SttG (Oij + flij + e^(2^ - Aj)x/-i(B)) • (2.14c) 

In ()2.14p we have used the indicator function X/-i(b) of the body to make 
explicit the support of p. Define T,Ti,Tij by 

Tf,^ = T{dt + ilJidx'f + 2Tjdx\dt + i^idx') + njdx'dx^ . (2.15) 

Then we have 
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Lemma 2.1. 

e''i2-^-phi,) = n,, (2.16a) 



e^^ = -r., (2.16b) 



e^{p + §^) = e-'^r + r/. (2.16c) 

For proof of Lemma l2 . 1 \ see appendix[Al After substituting (I2.16P into ()2.14p 
the reduced field equations take the forrrQ 

Af,U = 47rGx/-i(B)(e-^^T + r^^) - e^^'ojkiu;'' , (2.17a) 

D^e'^^Uij) = 87rGx/-i(B)Tj , (2.17b) 

Gij = 8'iTG{xf-i(B)Tij + Qij + Qij) . (2.17c) 

2.5. Integrability conditions. The quantities Qij and Qij satisfy the identi- 
ties 

SnGD^e^j = 2 {DiU)AhU , (2.18) 

and, using D^iUjj^ = 0, 

SirGD^nij = 2 [e^^ {DiU)uJkiuJ^^ - oJik D^^^^w/)] . (2.19) 
In case G 7^ 0, we obtain from ()2.17bp the integrabihty conditions 

D'Ti = Q, (2.20) 
together with the boundary condition 

Tin^l/-i{9B) = 0. (2.21) 

Further, we have 

D^nj - 2uj,jT^ = -{D,U){e-^^T + Tk^) , (2.22) 

and 

rijn^|/-i(ai3) = 0, (2.23) 

as a consequence of the contracted Bianchi identities for h.ij appUed to the left 
hand side of (I2.17cp . together with (f238]) . (f239]) and (I2.17ap . 

2.6. Implementing rotation. Define a vector field r]^ by 

r/^a^ = 7]'di . (2.24) 
The scalar product a = Qfiui^rf satisfies 

+ ^^7,^=0. (2.25) 



-2U 



a 



Since by assumption [J^,g) is stationary with respect to it holds that 77'' 
commutes with if and only if 77* does not depend on t. In particular, in the 



^Equation (|2.17bp corrects a typo in [31 (2.47)]. 
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case we are considering, the vector field rj'^ is itself a Killing vector if and only 
if the equations 

Cr^U = 0, (2.26a) 

Crji^i = , (2.26b) 

Crjhij = , (2.26c) 

hold. In these expressions the operator means the Lie derivative of the 
respective object with respect to r]^dk- Note that (|2.26b|) implies 

2ujijr]^ + A(e"^^a) = 0. (2.27) 

Define the velocity field by 

uf" = b-\e + ^V'') , (2.28) 

where is a real parameter corresponding to the rotation speed, and 6 is a 
normalizing factor, determined by n'^n^ = — 1 in f~^{B). It is important to 
note here that the rotational field + ^Irj^ in general will fail to be globally 
timelike for nonzero values of $7. However, for a suitable range of Q, it makes 
sense to require + Jlr/'^ to be timelike in the body. 

We now impose rotation of the body by requiring that the configuration 
satisfies the condition u^f^^/i = 0, i.e. 

/^,^(e^ + n7?'^) = 0. (2.29) 

Since T^yu'^ = —pu^, due to (|2.ip the stress energy tensor satisfies the relation 

n[^T,]X = 0. (2.30) 

It follows from (|230]) . that 

[uoTi^ - UiTof,)u^ = (2.31) 

holds, which, using (j2.29p . after some cancellations and multiplying by e~^^ 
gives 

(1 - Qe-^^afn + 0(1 - ne-'^^a)Tijri^ 

+ ne-'^^7]i[{l - ne-^^a)T + firj-ry^] = . (2.32) 

Equation (j2.32p can be proved by explicit computation, using (j2.16p and (j2.34p . 
As a consequence of (|2.32p we have 

Lemma 2.2. For sufficiently small Q, 

(1 - ne-^^a)Ti + ^Tijif + ne-^^r]i[T + f](l - ne'^^ a)-^Tj7]^ = . (2.33) 

2.7. Stress tensor. In order to write the field equations, we shall need the 
stress tensor for the elastic material. For consistency with the treatment of the 
static case considered in [T], we shall here make use of an analogous form of 
the stress tensor. Recall that assuming material frame indifference, the stored 
energy function e is a function of f^ and = f^^p.f^,v9^^- Taking equation 
(I2.29P into account, we find 

= /(^,,/^),,[e2^/i*^' + ine^^^p^rf + 02(e2^V''V'A: - e-^'')r)W] ■ (2-34) 
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In the computations below we shall make use of H"^^ defined by 7"^^ = 
f,2UjjAB^ so that 

Let 



Be 

^AB = -2-gj^, a^u = nf^,^f'',,aAB, = f"" ,^(TbcH^'^ ■ (2.35) 



It follows from the definition that 



(TAB = -2e 



-2U 



de 



AB ■ 



Our next task is to evaluate the dependence on ^1 of the terms occurring in 
the left hand side of (|2.16p . It is straightforward to verify that the following 
Lemma holds. 

Lemma 2.3. There are z,Zi,Zij depending smoothly an f^, gp,y, and their first 
derivatives, as well as 0, and G, such that the following equations are valid. 



e^(2^ - phij) = -e^{(Jij - ^Zij) , (2.36a) 



e^^ = -e^Oz., (2.36b) 



By the results of Lemma 12.31 and Lemma 12.11 we have 

Tij — 6 (<^2j ^^ij^ . 

We are now able to rewrite the integrability condition (j2.22|) in the form 

D^{e^aij) = e^(ne - a/)DiU + Q[D^{e^Zij) + 2e^ujijZ^ + zDiU] . 

Taking the above facts into account, we arrive at the system of equations 

AhU = 47TGxf-iiB)e^{ne - a/) + Qz) - e^Vz^*^' , (2.37a) 

D^(e^%j) = 8^Gx/-i(B)e^Jlzj , (2.37b) 

Gij = 87rG[-X/-i(B)e^((Tij - Qzij) + Qij + Qij] , (2.37c) 

D^{e^aij) = e^(ne - a/)DiU (2.37d) 

+ n[D^{e^Zij) + 2e^ujijZ^ + zDiU], mf-^{B), (2.37e) 

subject to the boundary condition 

{aij - nzij)n^\gj,_^^^^ =0. (2.37f) 
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2.8. Gauge reduction. Two of the equations in the system (|2.37p fail to be 
elhptic in the form given above, namely (I2.37b[) and (12.37cp . The reason for 
this failure is the related to the diffeomorphism invariance of the 4-dimensional 
Einstein equations. As in the static case, the method which shall be used to 
avoid this problem is make use of harmonic coordinates. 

Let □ denote the wave operator in {M,g). Taking into account of the fact 
that g^u is stationary, we have 

at = -^d^{g'"'^du)t (2.38) 

= e^^DVi • (2.39) 

Thus, e^^ D^ijji = precisely when the time t is harmonic. 
The left hand side of equation (j2.37b[) is of the form 

D'ie^'^uji,) = e^^[4D'Uu;ij + ^{A^jj - i?/Vfc)] " ^e^^DjD^ . (2.40) 

The term DjD^ipi causes this expression to fail to be an elliptic in ipi. However, 
the following reduced form of equation ()2.37bp . 

which modifies the left hand side by a quantity that vanishes if the harmonic 
time condition is satisfied, is elliptic in ■^j. 

Similarly, (j2.37cp fails to be elliptic due to the covariance of Rij. Following 
[H section 3.1], let V = h^^{T''-f^-f'--^) where fj-;. are the Christoffel symbols of 
a fixed Euclidean background metric on M. Then = is the condition for 
harmonic coordinates in M. By replacing Rij by Rij — D^^Vj^ we arrive, after 
rewriting equation (j2.37cp making use of the identity (3.11)] at the reduced 
Einstein equation 

- ^Ahhij + Qij{h, dh) = -SirGe^ {aij - hij a/ + Q{zij - hijZi''))xf-i(B) 

+ 2DiUDjU + e^^lhijUJMUj''^ - 2wifcw/] . (2.42) 

As in pp, we shall first solve the reduced system involving (|2.42p and (j2.4ip 
and once the solution is in hand show that the solution to the reduced system 
is actually a solution to the full system. We construct solutions by an implicit 
function theorem argument applied to a projected version of the field equations 
in material form. 

2.9. Field equations in material form. In the Eulerian picture, the do- 
main /~^(;B) depends on the unknown configuration /. This introduces a "free 
boundary" aspect in the Eulerian version of the field equations, which we will 
avoid by passing to the material, or Lagrangian form of the field equations. In 
this form of the equations, the configuration / is replaced by the deformation 
(j), and the entire system of field equations is moved to the extended body Mg. 
In particular, in this formulation, the elastic field equation lives on the fixed 
domain B. 
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The Piola transform of cij-^ is 

Similarly, we introduce the Piola transform of Zij. Since B has a smooth bound- 
ary, there is a linear extension operator which takes functions on B to functions 
on M|. In particular this allows us to define an extension cj) of cp which is equal 
to i outside a compact set. We use (j) to move the fields from space to M|, 
and use the bar notation introduced in [H section 3.2] to denote the quantities 
transported under (p. In particular, we define 



U = U o (j) , diU = diU o , = o (p ^ ^ij = hij o (f) . 

Note that for the barred quantities, it is the frame components which are pulled 
back, and not the tensor itself. Equation (j2.37ap in the material frame becomes 

'AiJi = AnGxBe^ine - a/) + ^z) - e^^Ui^ , in ]R| . (2.43) 
We remark that covariance of the Laplacian gives 

Next, equation (j2.4ip in the material frame becomes 

1- 



D'{e^u^..) + -Dj{e^Uj:).^i^ = 8nGxBe^nzj . (2.44) 



Equation (I2.42P becomes 



-Ahhij + Qij{h, dh) = -SirGe'^ {uij - hij ai + Q{zij - hijZi ))xb 



+ 2DiUDjU + e'^^lhijUJkiuo'''' - 2 - iOikUJjf"] . (2.45) 
Equations (|2.37ep and ()2.37fp become in the material frame 

n 



+ niDAie^Zi"^) + 2e^^^^ + zdiU] , inB, (2.46a) 

n 

subject to the boundary condition 

{a^^-nzi^)nA\gj^ = 0. (2.46b) 

2.10. Constitutive conditions. Similarly to the static case, we shall assume 
the existence of a relaxed reference configuration for the elastic material, which 
is such that suitable ellipticity properties hold for the elasticity operator evalu- 
ated in the relaxed state. The relaxed state is given by the body B, a compact, 
connected domain B C M| with smooth boundary dB, together with a refer- 
ence configuration i : W^. We assume a reference Euclidean metric 6 on 
M = is given. The body metric on is defined by = i*6. The 
relaxed nature of the reference configuration is expressed by the condition 



= 0, inB. 

(U=0,H=5b) 
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(2.48) 

{U=0,H=5b) 



The specific rest mass, i.e. the rest mass term in the relativistic stored energy 
function, should obey 

e(^) = e|(c/=o,H=<5e) - 
for some constant C > 0. Further, we assume that the elastic material is such 
that there is a constant C > such that the pointwise stability condition 

LabcdN^^'N^'' > C {6caSbd + 6cb6ad)N^'' N"^"" , in ^ , (2.47) 
holds, where 

LabCd{X) := y-QjjABQHCD ) 

In the isotropic case considered in this paper, e depends only on the invariants 
of 7"^^ = e^^ H^^ , defined with respect to the body metric {Sb)ab, cf. section 
12.11 It follows that e is independent of X and there are constants A and so 
that 

Labcd = ^Sab^cd + '^i^^C{A^B)D , (2.49) 
in terms of which the condition (j2.47p holds exactly when 

/i>0, 3A + 2/i>0, (2.50) 

cf. section 4.3]. The constants A and fi are apart from a common constant 
factor the classical Lame moduli. The inequalities (|2.50p are usually expressed 

by saying that the Poisson ratio defined by = 2(^') ^^^i^fy — 1 < < ^. In 

fact for most materials occuring in practice there holds j < u < |. 

We shall assume that the body is axisymmetric. To make this notion con- 
crete, let X* and X"^ be coordinates on and Mg, respectively, so that the 
Euclidean metrics 6 and 6b have components 6ij and 5ab, respectively. The 
body B is axially symmetric if there is a one-parameter subgroup of Euclidean 
motions, defined with respect to 5ab, which leaves B invariant. We may with- 
out loss of generality assume that the subgroup leaving B invariant is generated 
by the Killing field 

Tj'^dA = x'^di - x^d2 , (2.51) 

which necessarily is such that rj^ is tangent to dB. Given the axial Killing field 
r/^ on Mg, define a vector field r/* on M.^ by 

rj'di = U7]^dA) . (2.52) 

In particular, we may without loss of generality assume rj^di to be of the form 
rj^di = x^di — x^d2- 

In addition to the above mentioned conditions, we shall in the following 
assume that the elastic material is isotropic, cf. section [2Tl Recall that if the 
elastic material is isotropic, then A and hence also the stored energy function 
e depends only on the invariants Aj of 7"^^, defined with respect to the body 
metric 5b, cf. section [2Tl Consequently, in view of the discussion above, see 
in particular section 12.71 the reduced energy density p = ne can be viewed as a 
function p = p{Xi). 

The invariants A^ are functions of the form Aj = Ai(/, 5/, C/, ■i/'j, hij \ "q'^Cl). In 
the present case, we are using a coordinate system on B in which the metric 
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(5g has constant components, and hence the Aj do not depend on / but only 
on its derivatives. We may therefore write p as a functional p = r/, $7], 

where the symbol g is used as shorthand for the gravitational variables C/, 'i/'jj hij 
parametrizing the spacetime metric g^j^. 

3. Analytical setting 

We now introduce the analytical setting which will be used to construct 
solutions to the field equations. Fix a weight 5 S (—1, —\)- Further, fix p > 3. 
The parameters 6,p will be used to determine the weighted Sobolev spaces 
which are used in the implicit function argument. 

The system of equations in material form has the unknowns C/, V'*, hij. Let 

Br = W^'P{B) X Wf" X Wf" X e]'^, 

where E^'^^ is the space of asymptotically Euclidean metrics introduced in [H 
section 2.3], and let 

= [L^'iB] X B'-'/P'^dB)] X x x L^,_, . 

Thus, Bi is a Banach manifold and B2 is a Banach space. 

The residuals of equations (|2.46aP with boundary condition (|2.46b|) . (j2.43p . 
()2.44p . (I2.45p . which depend on the Newton constant G and the rotation velocity 
Q as parameters define a map : x i?i — > i?2- Thus, J- has components 
{J^(p,J^u,J^^,Th) corresponding to the components of B2, given by 



where 



(3.1a) 



of// 



n 



■dB 



— n[DA[e'' z. 



+ 2e 



n 



+ zdiU] , 



dB ■ 



and 



Aftf/ - 47rGxBe^(ne - ct/) + ^z) + e^^Wfc/w'^' , 
1- 



SttG xbg ^Zj , 
^hhij + Qij(h, dh) + SnOe^ {aij - hij ai + Q.{z 



V 



2DiUDjU 



(3.1b) 
(3.1c) 

(3.1d) 
(3.1e) 

hijZi))XB 

(3.1f) 



We now have T = J^i{G, O), (0, U, ^, h)). Write a general element of Bi as Z. 
We will use the implicit function theorem to construct solutions to J-" = for 
G,Vt close to G ]R2_ 

An essential assumption which allows us to introduce a relaxed configuration 
is that there is a reference Euclidean metric 5 onM = R^, and a diffeomorphism 
i : M| — > M|. As discussed in section [27101 an Euclidean metric on Mg is defined 
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by Si3 = i*6. Recall that B is assumed to be a connected domain with smooth 
boundary. 

From the constitutive conditions, cf. section [2. 101 we have that 

Zq = (i, 0, 0, 6ij o i) 

is a solution to the equation J-'{0, Zq) = 0. In order to apply the implicit function 
theorem at {0,Zq) it is necessary that the Frechet derivative D2J^{0, Zq) is an 
isomorphism. We see that J-'{0, Z) is of the form 
r - - Pi- £ 



^0(0, 



e^(e-^)9,C/, 
n 



- A I 



J^uiO,Z) = AhU + e^^ljkiuJ 



,kl 



Z) 



--Ahhij + Qij{h,dh) - 2DiUDjU - e [hijuJki^'''- - 2uJikUJj^] 



It follows from the constitutive conditions stated in section [2 . 1 Ul that D^!Ffi,(0, Z) 
is elliptic. 



3.1. Projected system. An analysis along the lines of 
that D2-^(0, Zq) is of the form 







section 4.2] shows 



V 



A 












where the entries are evaluated at Z( 






-\^) 

Q. The diagonal entries are isomorphisms 
between the weighted spaces given in the definition of Bi and i?2, with the ex- 
ception of D^Tfj). As in the static case this has a nontrivial kernel and cokernel, 
see the discussion in [1, section 4]. The kernel and cokernel can be identified 
with the space of Killing fields on {B^5b). Therefore, in order to construct 
solutions, we will consider the projected system 

where P/j : i?2 — *■ B2 is a projection operator which is defined exactly along 
the lines of [H section 4]. In particular, is defined to act as the identity in 
the second to fourth components of -621 while in the first component of B2 it 
acts as the unique projection along the cokernel of D^T^p^O, Zq) onto the range 
of Dfj,T^{0, Zq), which leaves the boundary data in the first component of B2 
unchanged. We use the the label B to indicate that operates on fields on the 
body and the extended body. We shall later need to transport the projection 
operator to fields on M^. 

Let {hi,Ti) denote pairs of elements in W^'P{B) x W^-^/'P^'P{dB). The restric- 
tion of to the first component of B2, which we here denote by the same 
symbol, is defined by setting Pg(6j,rj) = (6'j,rj), satisfying 



en 



dB 



(3.2) 
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for all Killing fields Pairs {b'i,Ti) satisfying this condition are called equi- 
librated. As discussed in ^1,, section 4], the definition of Pg implies there is a 
unique rji of the form r/j = ai + (3ijX^ , for constants Ui, (3ij satisfying f3ij = —Pji, 
such that 

b'i = bi- r]iXB ■ 

We further restrict the domain of PgJ^ to eliminate the kernel of D^FqJ^. 
By assumption, cf. section [2. 101 B has an axis of symmetry, which without loss 
of generality can be identified with the X^-axis. Fix a point Xq on the axis 
of symmetry of B, i.e. Xq has coordinates (0,0, X^) for some X^. Recah that 
the kernel of D(f,J^ consists of the Killing fields of {B, Sq). A killing field in B is 
determined by specifying its value and antisymmetrized derivative at one point. 
Following the proof of [1, Proposition 4.3], defin^ X to be the submanifold of 
Bi such that 

(</.-i)(Xo) =0, 8.nd6^,6ciAdB]{^-iy{Xo) = 0. (3.3) 

and define Y to be the range of the projection operator Pg. An application of 
the implicit function theorem to the map 

FbJ" : X ^ Y 

now gives the following result, analogous to [1, Proposition 4.3]. 

Proposition 3.1. Let T : Bi ^ B2 he map defined by \3. 1\) and let Pg be 

defined as in [H section 4.3]. Then, for sufficiently small values of Newton's 
constant G and the rotation velocity there is a unique solution Z = Z{G, Q), 
where Z = {<j), U, ipi, hij), to the reduced, projected equation for a self-gravitating 
rotating elastic body given by 

¥bH{G,^),Z)=0, (3.4) 

which satisfies the condition i3. 3\) . In particular, for any e > 0, there are 
G > 0,il > 0, such that Z{G, il.) satisfies the inequality 

||(/)-i||iy2,p(B) + \\hij -6ij\\^2,p + ||C/||,^2,p + < e- (3-5) 

The proof of proposition 13.11 proceeds along exactly the same lines as the 
proof of pj, Proposition 4.3] and is left to the reader. 

4. Equilibration 

Arguing along the lines of [1, section 5], we have the following corollary to 
Proposition 13. 1[ 

Corollary 4.1. For any e > 0, there are G > 0,Q > such that the inequality 
\\(p - i||iy2,p(B) + \\hij - 6ij\\^2,p + ||C/||w2,p + M\yy2,p < e (4.1) 

odd 

holds. 



The discussion here corrects some typos in the proof of W, Proposition 4.3], in particular 
the antisymmetrization in (|3.3p corrects the corresponding expression in [I]. 
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4.1. Eulerian form of the projected equations. Let Pj be the Euler- 
ian form of the projection operator, defined as in [H section 5.1] by 

P^-i(e)(n-(6o/))=n(Pe6)o/. 

Moving to the Eulerian form of the projected system, we find that we have 
constructed, for small G,Q a solution (cj), U, il^i, hij) of the following set of pro- 
jected equations, which it is convenient to write in terms of the stress energy 



components r, Tj , Tij . 

Pj-i(e)(D^ri, - 2L0,jT^ + {D,U){e-^^T + n'')) =0, (4.2a) 

njn^\df-HB)=0, (4.2b) 

AhU = 4^Gx/-i(6)(e-'^r + Tk') - c^'^'umu;'' , (4.2c) 

D\e''^uJij) + ^Dj{e''^D%) = 87rGx/-i(B)Tj , (4.2d) 

Gi, - + h,jDiV' = 87rG(x/-i(B)r,, + 9,,- + O,,) . (4.2e) 



Let Y = {f ,U,il)i,hij) be the Eulerian form of the solution to the projected 
form of the material field equations constructed in section [3?T1 From proposition 
I3.H the solution is unique. For the purposes here, we shall need to make 
the uniqueness property somewhat more explicit. An analysis of the proof of 
proposition 13.11 proves the following corollary. 

Corollary 4.2. Let the body domain B with metric 6js be given, with the corre- 
sponding background metric 6 on M = M^, and fix a point Xq in B and a vector 
field T]^ on B. Then the Eulerian form Y = {f^, U, ip, h) of the solution to the 
reduced projected system for a self-gravitating, rotating, elastic body defines a 
function of the form 

Y = Y{G,n;[B,6Bj,Xo,ri]). 

4.2. Equivariance. We now analyze some of the consequences of the consti- 
tutive conditions imposed in section 12.101 Recall that in particular, in view 
of frame indifference and homogeneity, and the discussion in section 12. 10^ the 
reduced stored energy function p is of the form 

where f^ is the configuration, g is used as shorthand for the fields U, ipi, hij on 
M parametrizing the spacetime metric 5^;^, and r/* is the axial vector field on 
M specified in section 12.101 Let u be a spatial diffeomorphism, i.e. t o a = t. 
Then by frame indifference (i.e. general covariance) we have 

p[f o <^*g; c^*??, ^ = p[f, 9;v,^]°c^ ■ (4-3) 

Further, as a consequence of the isotropy of the elastic body, for any isometry 
E of {B,Sb), we have 

p[j:of,g;r],n]=p[f,g-r,,n]. (4.4) 
The transformation properties stated in (14. 3p and (14. 4p give the following Lemma. 
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Lemma 4.3. Let {f^,U,il)i,hij) be as in corollary Let H he a diffeomor- 
phism of B leaving the data {Xo,6ab,V^) invariant. Then the diffeomorphism 
a of M defined by requiring that ioT, = a o i on all of M| is an isometry in the 
sense that it leaves all of {f'^,U,il^i,hij,r]'') invariant. 

Proof. First note that a is by construction an isometry of the flat background 
metric 5 entering the projected, harmonicahy-reduced field equations and that 
(c*??)* = 1]^ trivially from the construction of ry*. Using these facts together 
with the equivariance property expressed in (j4.3j) and ()4.4p we have that 

((S-i o / o a)^, a*U, {a*h)ij) (4.5) 

is a solution with the same data. By the uniqueness property made explicit in 
corollary 14.21 we then have 

((E-i o / o a)^, a*U, {a*i^)i, {a*h)i,) = (/^, U, i^i, hi,) . (4.6) 

□ 

If a is as in lemma then we also have 

a*T = T , {cr*T)i = Ti , {a*T)ij = nj . (4.7) 

Lemma 14.31 has the following corollary which will play an important in the 
proof of orthogonal transitivity, see section 14.51 below. 

Corollary 4.4. Let {f"^, U,ipi, hij) he as in corollary \4.S\ Let S be an isometry 
of (Bjds) such that T,{Xq) = Xq and (Y^^rj)^ = —rj^, and let a be a diffeomor- 
phism of M such that ioT, = a o i on all o/Mg. Then a is an isometry of hij 
and we have 

(S-i of^oa, a*U, (a»„ {a*h),,) = (/^, U, -^Pi, h,,) . (4.8) 

Proof. The transformation ipi —tpi, rj^ — > —r]^ leaves H^^ and hence all the 
field equations invariant. Therefore it maps a solution to another solution. By 
uniqueness it follows that the solution with data B, Sq, 6, —rj^, G, is given by 
(/^, U, -ipi, hij). The result follows. □ 

Recall that the reference state is axially symmetric, i.e. rj^ is an axial Killing 
vector in Euclidean space leaving B invariant. Denoting by S the flow of ry"^ 
and correspondingly using a to denote the flow of rj'', we have the following 
infinitesimal version of of Lemma 14. 3[ 

Lemma 4.5. Assume that B is axially symmetric with axial Killing field rj^, 
as discussed in section [2.1(A Then 

/^(x)7?Xx)=r?^(/(x)), (4.9) 

%. c 

n'd, = f*{v^dA). (4.10) 

and 

CnU = 0, (4.11a) 

/:„V'i = 0, (4.11b) 

Crjhij = 0. (4.11c) 
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By the antisymmetry of ujij we have, using ()2.27p . that 

£^(e-2^a) = 0. (4.12) 
Furthermore, from (14.7jl apphed to the flow of rj^, we infer that 

£„r = 0, £„ri = 0, £^Ti, = . (4.13) 

4.3. Divergence identities. Now turn back to Eq. (j2.33p . Taking the diver- 
gence of this equation and using (j4.11al4.11cl4.12l4.13|) and (j4.2ap . gives 

= (1 - ne-^^a)DWi - nD'{e-^^a)Ti 

+ Or?^'(I/-i(B) - P^-i(e))[DV,, - 2u;,iT^ + {D,U){e-^^T + r^^^)] (4.14) 

use (|2.26bp 

= (1 - ne-^^a)DWi 

(4.15) 

It also follows directly from (I2.33P and the fact that rf is parallel to the boundary 
of f~^{B) that the boundary condition 

nn'\f~i(dB) = (4.16) 

holds. Let W = e'^^ D^ipi. The first term in the left hand side of ()4.2dp is the 
divergence of a 2-form, and therefore its divergence is zero. Hence, taking the 
divergence of both sides of ()4.2dp , and using the fact that (|4.16p holds for the 
case of an axisymmetric body, gives the identity 

AhW=167rGxf^iiB)D'n. (4.17) 

Equation (j4.15p gives the form of the right hand side in (j4.17p . Let 

LV, = AhVi + R^^Vk , (4.18) 

and note 

Using (j2.18p and ()2.19p we find after taking the divergence of both sides of 
(gUD, when G / 0, that 

LVi = -m7TGxf-i(B)[D'Ti, + (A^7)(e-^^r + r^^^)] + Aoj.kD^e^^ co/) 
use (l42d]) and i^Alai\ 

= -levrGpV,, + {DiU){e-*^T + Tk'')]Xf-HB) + M,[8TrGxf-iiB)r' - \d^W] 
+ lQTiG¥f-,^js)[D'Tij + {D,U){e-^^T + r^'^) - 2a;ijr^']x/-i(B) 

= -16^G(I^-i(e) - F^-i(B))[DV,, + (AC/)(e-^^T + Tfc'^) - 2u:i^T^]xf-^B) 
- 2uJijD^W . 
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Let 

= -16TTG[D^n, + {DiU){e-^^T + r,.^) - 2u;ijr^] . 
Then we have the following system of equations for W, Vi , 

(1 - f7e-2^a)AT^ = - F/-i(e))2,X/-i(e) , 

LV, = - P;-i(e))^^X/-i(B) - 2u:i,D^W. 

Arguing as in the proof of [1, Lemma 5.7], we have that 

(Ej-i(e) - Pj-i(e))2:iX/-i(B) = n(Ci o /)x/-i{B) , 
for some Q which is a Killing field in B. Hence we have the equations 

(1 - ne-^'^a)^W = ^rfn{Cj o f )xf-^B) , (4.19a) 
LV^ = n{Ci o /)x/-i(B) - 2u;ijDm . (4.19b) 

4.4. Main theorem. We are now able to prove the following 

Theorem 4.6. For sufficiently small values of G, il., with G non-zero, the so- 
lution to the reduced, projected system of equations for a stationary, rotating, 
elastic, self- gravitating body (3^, is a solution to the full system of equations 
{2.11^ for a stationary, rotating elastic, self- gravitating body, together with the 
integrability conditions of section \2.5[ In particular, this solution corresponds 
to a pair {f^,gp.v), which solves the full Einstein equations G^^ = SirGT^^^. 

Proof. Using the estimate of Corollary 14.11 and the multiplication properties of 
the weighted Sobolev spaces, cf. [H section 2.3], one checks that 

u^^l(^wl\, n^,ew^f_,, e,,ew^f_,, 

with corresponding estimates. Hence we find, using equation (14.2eP for hij, in 
the equivalent form (12.42p . equation (14.2d|) for ipi, making use of equation (j2.40p 
to express it in a form suitable for estimates, as well as equation ()4.2cP for U, 
that the conclusion of [1] Lemma 5.2] for hij holds also in the present case, 
namely 

hij = 5ij + — + /l(2) ij ) 

for constants 7jj, with /i(2)ij £ ^2&' ■ sufficiently small G,VL we have the 
estimate 

11^(2) ijWw^f + IItII - C'dl^ii - ^ijWwf^ + 11*^ ~ i|lvi/2,P(B)) • 

For brevity, we shall in the following write estimates of the above form using 
\\Z — ZqWbi where the norm refers to that induced from the Banach spaces 
using in defining the space Bi , cf . section [3l We shall further write inequalities 
of the form a < Gb where C is a constant which is uniformly bounded for small 
G,fl as a < b. 

Given this result about the asymptotics of hij, the conclusion of [H Lemma 
5.4] concerning Vi holds, and hence also the partial integration result [U Lemma 
5.5] and the estimate of Lemma 5.6]. Now define the operator Q : L^_3(]R^) — j 



22 L. ANDERSSON, R. BEIG, AND B. SCHMIDT 

as in [T, section 5.2]. Given a basis {C(k)}k=i for the space of Killing fields, 
we set 

Q>iizi) = D^idf^h, K = l,...,6. 
M 

Since W = e^'^D^ipi, the term uJ^jDm in (l4T9bl) satisfies UijDm G L^^.g 
and we have the estimate 

||u;.,Z)^'H^||iP_^(K3) < ||Z-Zo|bJ|n(Co/)||^p(^.,(^)). (4.20) 

Recall that from the construction of Q we have for small G, il., the equivalence 
of norms 

I |Qn(C o /)x/-i(B) I |r6 < I Id |iR6 < I |Qn(C o /)x/-i(B) I Ir6 , (4.21) 
where if = a* + P^jX^ , IICIIrs is defined by 

\\C\\le=Y.iar + Y.(f^',?. 

i i<j 

Due to the properties of Q, the analogue of (j4.2ip holds also for | |?i-(C°/)l 
Applying Q to both sides of ()4.19bp . we have using (jOO]) and (fOT|) . 

I|C||r6<||QLF||k6 + ||Q^W||m6 

<||QLF||K6 + ||Z-Zo|bJ|C||R6, 

and hence 

||C||R6<||QLy|b6. (4.22) 

Recall that for G, Q sufficiently small, we also have due to Corollary 14.11 that 
\\Z — ZqWbi small. We now have the chain of inequalities for G,Q sufficiently 
small. 



HC ° f)\\LP{f-HB)) + W^DWWlp 



' S-1 

< 

<\\C\y + \\z-Zo\\B, 
<IICIIr6 

use K22\i 

<\mv\y. 

By the inequality proved in [1, Proposition 5.8] we have 

WQLVW^. <\\Z - Zo\\bA\V\\^2^,^, 
which together with the above gives 

\\y\\w?- <\\Z-Zo\\Bmw^- ■ (4-23) 

o — 1 o — 1 

By choosing G,Q sufficiently small, we can make \ \Z — Zo\\bi small enough so 
that (I4.23P gives the inequality 



l^lll^/2,p < 



1 
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which imphes 



V = 0. 



Due to the vanishing of V, it follows from ()4.22p that also C = 0; and hence we 



This means that the solution of the projected system of equations (14. 2p is actu- 
ally a solution to the full system of field equations (I2.17P for the rotating elastic 
body, together with the integrability conditions discussed in section 12.51 

It remains to demonstrate that the solution {f^,U,il>i,hij) to (I2.17P con- 
structed in this proof corresponds to a Lorentzian spacetime {M , g^y ) solving 
the Einstein equations for the elastic body. The solution we have found yields 
via (12. 4p a Lorentz metric g^i, at some time to together with its vanishing first 
and second time derivatives at Iq, as well as a configuration together with 
its non- vanishing first time derivative at Iq. These solve the Einstein equations 
at to- We extend the spacetime metric off to by requiring it to be t- independent 
and by requiring it to satisfy (|2.29p for all times. This constructs a space- 
time {M,g^u) which is axisymmetric and stationary and a configuration which 
is axially symmetric and helical. Thus, by the discussion in section [T2t the asso- 
ciated energy momentum tensor is time independent. This shows that (^J^,g^i/^ 
together with the configuration provide a solution to the full Einstein equa- 
tions. □ 

We remark that the solutions we have found are static exactly when Q = 0. 

4.5. Orthogonal Transitivity. Let {M,g^y) be a stationary spacetime con- 
taining a rotating elastic body as constructed in Theorem 14.61 We have shown 
in section |4] that (Ai,g^y) admits a two-parameter, abelian group of isometries, 
generated by the Killing fields ^,^,7]^. In fact, since r]^^ is the pullback of the 
axial vector field acting on the body, the group can be taken to be the cylinder 
R X S*^ . The question arises if this group acts orthogonally transitively on , as 
is the case for perfect fluids. Recall that a group acts orthogonally transitively 
if the the distribution perpendicular to the generators of the group action is 
Frobenius integrable. 

Define w^^x = S^ipVuCx] and let lu'^j^^x be defined with respect to r/^ in the 
analogous manner. Orthogonal transitivity is equivalent to the conditions 



see [HI (2.53)]. The spacetimes constructed in this paper have metrics which 
fail to be smooth at the boundary of the body f~^{B). 



have 



W = 0. 



(4.24a) 
(4.24b) 



Proposition 4.7. Let {M,g^y) he a stationary spacetime containing a rotating 
elastic body as in Theorem \4.6[ , with stationary and axial Killing fields ^^,r]^- 
Then, if Q > is sufficiently small, equation \4-24^ holds in M.. 
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Proof. The conditions (j4.24p can be restated in the space manifold M as 

e^^7][,u;jk]=0, (4.25a) 

e~'^^r][iDjr]k] + a^iUJjk] = . (4.25b) 

Here LOij = d^iipj] and a = C'^r]^, as above. Equations [3l (2. 60), (2. 61)] are 
equivalent to (|4.25|) but written in terms of a different representation of the 
spacetime metric. 

By [31 equations (2.51)-(2.52)], we have 

4V'^(r?[^u;^,A]) = (^eR^.[uCxr]p] , (4.26a) 

4V'^(e[p'^V-A]) = -Qv^^R^luCxrip] . (4.26b) 

It follows that 

^^ivip^pux] - ^^[p^'pux]) = l27rG{e' + ^t)T^wixr^p] ■ (4.27) 

Since the velocity vector = b~^{^'^+^r]'^), see equation (I2.28p . is an eigenvec- 
tor of the stress energy tensor, cf. equation ()2.ip . the right hand side of ()4.27p 
is zero. The left hand side of (I4.27P is the divergence of a 4-form, i.e. in terms 
of the exterior derivative and the Hodge dual, we have an equation of the form 
•k d-ka = 0. In particular, -ka is a scalar function which is constant, d*a = 0. 
In the situation under consideration, a vanishes on the axis = = 0, and 
hence it is zero everywhere. In this argument we made use of the fact that 
is well defined in all of f~^{B). This holds for sufficiently small values of Q, 
since then the vector field + Qr]'^ is timelike in all of f~^{B). In terms of the 
space manifold M, we have shown that 

- e^^{l - Qe-^^a^iLo^k] = ■ (4.28) 

This can of course also be checked directly from the three dimensional field 
equations. Note that relation ()4.28p becomes vacuous in the static case = 0. 

By the above argument we have shown that the two equations ()4.25p are 
linearly dependent if 7^ 0. Thus, in order to show that both equations in 
(j4.25p hold, it is sufficient to show that r][iDjr]j^j = 0. To see this we argue as 
follows. It follows from the axisymmetry of the body that there is a discrete 
isometry S of {B,6b), consisting of reflections in planes containing the axis, 
which maps rj^ to —rj^. An explicit choice of S is given by 

By corollary 14.41 the construction of r/* , we have that the diffeomorphism 
o" of M defined by S o i = cr o i is an isometry of hij, which has the property 
that {(T*riY = —rf. We can now conclude that reflections at planes through the 
x^-axis preserve both U and hij and send both ipi and "rf to their respective 
negatives. So in particular these transformations preserve vectors tangent to 
these planes, and since they send rf to —rf" and preserve inner products, rf has 
to be orthogonal to these planes. Consequently rf is hypersurface orthogonal, 
i.e. 

r?[iD,%]=0. (4.29) 
It follows, using (1428]) . that holds. □ 
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Remark 4.1. Recall the identity valid for Killing vectors 

3D\rj^Dj7],]) = 2rjyRk]iv'- (4-30) 
Inserting !i4.30\ ) into \2.14c^ , using \2.27^2.26a\) and finally U.28^ , there results 

2,&[{l - fie-2^a)i77[il),%]] = 167rG(l - ^e'^^ a)\r,%^rik^ . (4.31) 

Thus we have inferred that rf is an eigenvector of the stress tensor. This latter 
fact could have also been shown directly from the reflection symmetry without 
using the Einstein equations. 

Remark 4.2. In the case of a smooth spacetime, it follows from 1^4 ■24\ ) o.'f^d the 
Frobenius theorem that the distribution perpendicular to ^'^,7^^ is integrable, in 
the sense that there are smooth 2-surfaces inM orthogonal to the span of^^^^rj^. 
The spacetimes constructed in Theorem \4.6\ have in this paper shown to be Wf^^. 
Although the spacetimes containing a rotating body can in fact be shown to be 
real analytic away from the boundary of the body, f^^{dB), a further analysis 
is needed to show that an appropriate version of the Frobenius theorem applies. 
This question will be studied in a later paper. 



Appendix A. Proof of Lemma 12.11 
We have A = e^^ p. From 



- 



we get 

Using the form of g^^ , cf. (|2.5p . we have 



dg- 



dhv 



Define T,Ti,Tij by 

T^^ = T{dt + ilJidx'f + 2Tjdx\dt + il^idx") + njdx'dx^ 

Then, 

Tij = Tij + 2^i7pj-j + Tlpiljjj , 

Toi =Ti + T1pi , 

^00=^", 
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We calculate 

u^2^ - phi,) = e-2^(2^ - KK,, 



Tij - Tooipiipj + A{gij - gooipitpj) - Ae~^^/iy 



dip'^ dtp 

-2U 



e 



OA dg^''' 



dgi^" dtjj 



= -Tio + ipiToo - Agio + A^j^oo 

= -n, 

= e-'^{l{T,, + Ag,,)^-2A) 

= e-2^((T^, + A5;..)(5^^ + 2e-^''S%d%) - 2A) 
= e-2^(r/ + 4A + 2e-2^Too + le-'''^ Ag^o - 2A) 
= e-2^(r/ + 2e-2^Too) 
= e-^^r + r/ . 
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